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Pore-scale  simulations  of  flow  and  transport  through  a  filter  were  used  to  predict  clean-bed  filtration  of 
nanoparticles  from  a  3D  image  of  the  filter.  X-ray  micro-computed  tomography  was  used  to  obtain  the 
geometry  and  topology  of  the  filter  consisting  of  flattened  half-spherical  collectors.  A  Lattice-Boltzmann 
method  was  used  to  model  fluid  flow  and  particle  transport  in  the  volumetric  image  of  the  filter.  Nanopar¬ 
ticles  are  considered  that  are  so  small  that  diffusion,  as  compared  to  the  processes  of  gravitational 
settlement  and  interception,  dominates  particle  deposition.  A  correlation  equation  is  determined  for 
the  average  collector  efficiency  for  diffusion  through  regression  analysis  performed  for  a  set  of  numerical 
breakthrough  experiments  for  a  range  of  suspended  particle  sizes  and  Darcy  velocities.  This  correla¬ 
tion  quantifies  the  collector  size  by  the  surface  average  equivalent  sphere  diameter.  The  newly  derived 
correlation  agrees  well  with  experimental  data. 
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1.  Introduction 

Nano-sized  particle  transport  provides  solutions,  concerns 
and  potential  hazards.  For  instance,  the  use  of  manufactured 
nanoparticles  in  groundwater  remediation  is  a  new  and  promising 
technology  1 15,47].  The  discharge  of  nanoparticles  during  reme¬ 
diation,  however,  poses  a  potential  risk  to  public  health  [25,36 1. 
Moreover,  the  transport  of  viruses/pathogens  with  sizes  on  the 
order  of  nanometers  in  saturated  porous  media  is  receiving  more 
and  more  attention  [34,3,37].  These  facts  demonstrate  the  need 
to  investigate  the  fate  and  transport  of  nanoparticles,  including 
particle  deposition,  in  natural  and  engineered  porous  media. 

Particle  deposition  has  been  conceptualized  to  occur  in  two 
steps  [13,12]:  (1)  transport  of  the  suspended  particles  to  a  collec¬ 
tor  surface  by  advection  and  diffusion  as  quantified  by  the  collector 
efficiency  rjo  [-];  and  (2)  attachment  of  the  particles  to  the  collector 
surface  by  interaction  forces  as  quantified  by  the  collision  efficiency 
a  [-].  One  can  use  a  and  r) 0  to  model  the  deposition  term  in  the 
following  macro-scale  equation  for  particle  transport  in  a  sphere 
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packing  during  clean-bed  filtration  [17]: 

^  +V  .(Dc)  =  V.(DmVc)-  (1) 

where  c  [kg/m3]  is  the  particle  concentration,  t  is  the  time,  v  [m/s] 
is  the  pore-water  velocity,  Dm  [m2/s]  is  a  tensor  that  accounts 
for  hydrodynamic  dispersion  and  molecular  diffusion,  e  [-]  is  the 
porosity,  and  dc  [m]  is  the  grain  (collector)  diameter.  The  first  term 
on  the  left  hand  side  of  the  equation  describes  particle  accumula¬ 
tion.  The  second  term  describes  advective  transport.  The  first  term 
on  the  right  hand  side  of  the  equation  describes  particle  transport 
that  does  not  occur  with  the  mean  pore-water  velocity,  and  the 
second  term  describes  particle  depositioa 

Because  of  the  computational  cost  of  pore-scale  simulations  of 
particle  transport  in  a  representative  elementary  volume  (REV)  of 
a  filter,  unit-cell  approaches  have  been  traditionally  employed  to 
derive  correlation  equations  for  rf0.  Such  unit-cell  approaches,  for 
example,  represent  a  filter  of  spherical  collector  particles  by  (1 )  an 
isolated  spherical  collector  that  is  surrounded  by  a  laminar  flow 
field  of  infinite  extent  [45],  or  (2)  Happel’s  sphere-in-cell  model 
[  1 9],  a  single  collector  that  is  surrounded  by  an  imaginary  spherical 
shell  in  which  laminar  flow  takes  place  with  hypothetical  boundary 
conditions  for  fluid  velocity  and  particle  concentration  [33,40]. 

The  aforementioned  unit-cell  approaches  have  the  following 
shortcomings: 
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Nomenclature 

List  o/symbo/s 

c  particle  concentration  (kg/m3) 

Cn  distribution  function  for  particle  concentration 
dc  diameter  of  collectors  (m) 

dp  diameter  of  particles  (m) 

Dc  surface  average  equivalent  sphere  diameter  (m) 

Dm  diffusion/dispersion  tensor  (m2/s) 

Dp  particle  diffusion  tensor  (m2/s) 

/n  distribution  function  for  fluid  density 

g  gravitational  acceleration  (m/s2) 

k  Boltzmann  constant  (J/K) 

L  height/length  of  the  column  (m) 

NPe  Peclet  number 

Nr  interception  number 

p  pressure  (Pa) 

Sc  surface  area  of  a  collector  (m2 ) 

At  time  step  (s) 

T  temperature  (K) 

U  superficial  velocity  or  Darcy  velocity  (m/s) 

v  pore-water  velocity  (m/s) 

t'o  Ax/At  (m/s) 

Vc  collector  volume  (m3 ) 

Ax  lattice  spacing  (m) 

a  collision  efficiency 

rjo  collector  efficiency 

r)D  diffusion  collector  efficiency 

rjcon  collector  efficiency  from  correlation  ( 1 6) 

rjsim  collector  efficiency  from  numerical  simulations 

r)eXp  collector  efficiency  from  experiments 

e  porosity 

v  kinematic  viscosity  (m2/s) 

a  surface  charge  density  (fxC/cm2) 

pp  particle  density  (kg/m3 ) 

pf  fluid  density  (kg/m3) 

tf,  x c  relaxation  times  for  fluid  and  particles,  respectively 


1.  Correlation  equations  based  on  single-collector  approaches  do 
not  accurately  account  for  the  complex  flow  and  concentration 
fields  in  real  porous  media,  because  real  pore  geometries  and 
topologies  are  much  more  complicated  than  those  of  unit  cells. 

2.  Even  if  a  correlation  equation  based  on  a  single-collector  model 
is  deemed  to  be  sufficient,  existing  single-collector  models  are 
limited  in  their  ability  to  describe  filtration  in  packings  of  non- 
spherical  collectors,  with  a  notable  exception  in  which  [351  used 
a  unit-cell  approach  to  study  the  influence  of  collector  shape  and 
roughness  on  p0* 

In  order  to  better  represent  the  topology  of  a  filter  pore  space, 
simulations  in  pore  networks  have  been  used  to  investigate  filtra¬ 
tion  in  porous  media  [6,38].  These  simulations  partition  the  pore 
space  into  smaller  elements  (pore  bodies  and  pore  throats)  and 
simulate  transport  in  each  of  these  elements  using  simple  alge¬ 
braic  rules  for  the  incoming  and  outgoing  particle  concentrations. 


It  is,  however,  difficult  to  make  quantitative  predictions  from  a  3D 
image  of  a  porous  medium  using  current  pore-network  modeling 
approaches  for  particle  filtration,  because  they  do  not  represent  the 
topology  of  the  porous  medium. 

Due  to  continuously  increasing  CPU  speeds  and  the  availabil¬ 
ity  of  massive  parallel  computers,  Lattice-Boltzmann  (LB)  methods 
[43,10]  have  become  attractive  for  performing  pore-scale  simula¬ 
tions  of  a  broad  range  of  transport  phenomena  in  porous  media. 
The  LB  method  addresses  complex  pore  geometries  within  parallel 
computational  platforms,  which  in  turn  allows  the  simulation  of 
large  domains.  LB  methods  have  been  used  to  simulate  fluid  flow 
(e.g.  [30[)  and  small  particle  transport  (e.g.  [23])  in  porous  media. 

Few  LB  studies  have,  however,  investigated  filtration.  Przekop 
et  al.  [32]  studied  particle  deposition  in  fibrous  filters.  Lin  and 
Miller  [26]  investigated  the  flow  field  during  filtration.  Baumann 
and  Werth  [4]  compared  flow  paths  and  particle  velocities  in  a 
micro-model  to  a  2D  LB  model.  Chen  et  al.  [7]  studied  the  per¬ 
meability  reduction  due  to  colloid  deposition  by  combining  X-ray 
tomography  and  pore-scale  LB  simulations.  Long  and  Hilpert  [27] 
developed  a  correlation  equation  for  an  average  single-collector 
efficiency  for  clean-bed  filtration  in  sphere  packings. 

The  objective  of  this  paper  is  to  demonstrate  that  pore-scale 
simulations  (and  in  particular  an  LB  method)  of  flow  and  transport 
can  be  used  to  predict  clean-bed  filtration  of  nanoparticles  from 
a  3D  image  of  the  filter,  which  consists  of  collectors  that  are  too 
non-spherical  to  apply  correlation  equations  for  rj0  that  have  been 
developed  for  spherical  collectors  (e.g.  [33,40,27]).  For  instance, 
clay  particles  are  not  spherical,  and  sand  and  gravel  particles  are 
not  perfect  spheres.  Also  manufactured  filters  do  not  always  consist 
of  spherical  collector  particles.  This  work  considers  collectors  that 
resemble  flattened  half  spheres.  We  will  further  consider  filtration 
of  nanoparticles  that  are  so  small  that  diffusion,  as  compared  to  the 
processes  of  gravitational  settlement  and  interception,  dominates 
particle  deposition.  Hence,  tj0  can  be  approximated  by  the  diffusion 
col  lector  efficiency  rfD.  In  the  experimental  portion  of  this  paper,  (1) 
a  column  filled  with  flattened  half  spheres  is  imaged  using  X-ray 
micro-computed  tomography  to  obtain  the  geometry  and  topol¬ 
ogy  of  the  packing  in  3D,  and  (2)  clean-bed  filtration  experiments 
are  performed  in  this  column  using  nanoparticles.  In  the  theoreti¬ 
cal  portion  of  the  paper,  (1 )  a  LB  method  is  used  to  model  particle 
transport  in  the  3D  image  of  the  filter,  (2)  a  correlation  equation  for 
r}D  for  nanoparticle  filtration  is  developed  from  a  regression  anal¬ 
ysis  of  a  set  of  numerical  breakthrough  experiments,  and  (3)  the 
prediction  of  the  correlation  equation  is  compared  to  the  column 
experiments  obtained  for  a*  1.  We  also  comment  on  the  column 
experiments  that  were  performed  under  unfavorable  attachment 
conditions,  a  <  1. 

2.  Experimental  methods 

2. 1 .  Particle  suspension 

For  the  nanoparticles/suspensions,  we  used  four  types  of 
surfactant-free  polystyrene  latex  particles  (Interfacial  Dynamics 
Corporation,  Portland,  OR),  each  of  which  has  a  relatively  uniform 
size  distribution  (see  Table  1).  The  120-nm  Latex  A  particles  con¬ 
tain  amidine  groups  on  their  surfaces,  and  therefore,  have  a  positive 


Table  1 

Some  properties  of  the  nanoparticles  used  in  our  experiments. 

Label  Diameter  from  TEM  by  manufacturer  Diameter  from  PCS  from  our  measurements  Surface  charge  density,  a 

Latex  A  I20±l7nm  134nm  +8.1  p,C/cm2 

Latex  B  19  ±  3  nm  28  nm  -2.4piC/cm2 

Latex  C  43  ±  7  nm  57  nm  -1.7^C/cm2 

Latex  D  93±5nm  109nm  -2.2^C/cm2 
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surface  charge  density  or.  The  Latex  B,  C,  and  D  nanoparticles  con¬ 
tain  sulfate  groups  and  have  negative  surface  charge  densities  in 
the  entire  pH  range  used  in  this  study.  The  particle  density  pp 
is  1.055  g/cm3  for  all  particles  as  reported  by  the  manufacturer. 
Particle  diameters  were  measured  by  the  manufacturer  using  trans¬ 
mission  electron  microscopy  (TEM).  These  values,  however,  are 
smaller  than  the  ones  obtained  in  this  study,  which  uses  photon 
correlation  spectroscopy  (PCS)  of  suspensions  of  the  nanoparticles 
in  distilled  water  (Brookhaven  Instruments  Corporation).  In  PCS, 
the  hydrodynamic  diameter  dp  is  inferred  from  the  measurement 
of  the  mean  diffusion  coefficient  of  the  suspended  particles  and  the 
Einstein-Stokes  equation.  We  used  the  particle  sizes  determined  by 
PCS  as  model  inputs  due  to  their  relevance  to  particle  transport  in 
aqueous  phase. 

2.2.  Porous  media 

About  uniform-sized  non-spherical  beads  with  an  approxi¬ 
mate  shape  of  flattened  half  spheres  (supplied  by  Atofina/Arkema 
chemical  company)  are  used  as  the  collectors  in  the  deposition 
experiments.  The  about  half-sphere  shape  is  due  to  the  manufactur¬ 
ing  process,  during  which  soft  spherical  PVDF  beads,  created  from 
a  PVDF  melt,  deformed  and  flattened  due  to  the  action  of  gravity 
before  complete  solidification.  The  beads  are  made  of  polyvinyli- 
dene  fluoride  (PVDF,  Kynar@  460).  Polyvinylidene  fluoride  resin 
was  chosen  due  to  its  relevance  in  environmental  technology  such 
as  membrane-based  water  treatment.  The  density  of  the  PVDF 
beads  measured  by  the  manufacturer  is  1.688  g/cm3.  The  size  in 
the  largest  dimension  is  about  5  mm.  The  surface  charge  of  the 
PVDF  beads  was  not  determined  in  this  study  but  is  expected  to  be 
slightly  negative  according  to  streaming  potential  measurements 
on  a  microfiltration  membrane,  which  was  fabricated  from  simi¬ 
lar  materials,  in  distilled  water  [24].  A  X-ray  photon  spectrometry 
(XPS)  analysis  was  conducted  and  showed  that  the  surface  com¬ 
position  of  Kynar®  461  (powder  form  of  Kynar@  460)  is  almost 
exclusively  C2H2F2.  Oxygen  is  not  detected  at  a  noticeable  level, 
indicating  the  absence  of  hydroxyl  or  other  oxygen-containing  ion- 
lzable  groups  on  the  surfaces  of  the  beads.  The  beads  were  rinsed 
extensively  with  ultrapure  water  generated  by  a  MilliQ  water 
system  (Millipore)  before  the  filtration  experiments.  In  order  to 
explore  favorable  chemical  conditions  (a  =  1 )  for  Latex  A  particles, 
deposition  experiments  were  also  performed  using  PVDF  beads 
which  were  immersed  in  a  solution  with  sodium  dodecyl  sulfate 
(SDS)  for  4h,  and  the  results  were  compared  to  experiments  using 
PVDF  beads  without  SDS  coating. 

23.  Solution  chemistry 

In  all  of  the  experiments,  the  buffer  solution  consists  of  1  mM  of 
CaCl2  to  control  ionic  strength,  and  0.1  mM  of  NaHCCh  to  maintain 
the  pH  at  7.5.  For  the  tracer  breakthrough  experiments,  NaN03  is 
added  at  a  concentration  of  0.1  mM  as  a  tracer.  All  chemicals  used 
are  reagent  grade. 

2.4.  Deposition  experiment 

A  cylindrical  column  with  an  inside  diameter  of  2.54cm  and  a 
height  of  29.3  cm  is  used  in  the  deposition  experiments.  Fresh  PVDF 
beads  are  densely  packed  in  this  column  prior  to  each  experiment. 
The  measured  average  porosity  of  the  entire  bead  packing  was  43%. 
Particle  suspensions  with  a  concentration  of  2  mg/L  were  prepared 
by  adding  latex  stock  suspensions  to  the  aforementioned  buffer 
solutions  and  sonicating  them  for  1  h  to  break  up  any  aggregates. 
Aggregation  of  latex  particles  at  different  calcium  concentrations 
was  studied  in  a  series  of  jar  tests.  Addition  of  1  mM  of  CaCl2  did 
not  cause  noticeable  turbidity  increase  in  the  time  frame  of  the  fil¬ 


tration  experiment,  indicating  stability  of  the  particles  with  respect 
to  aggregation.  The  packed  column  was  first  flushed  with  five  pore 
volumes  of  ultrapure  water,  followed  by  flushing  with  another  five 
pore  volumes  of  the  buffer  solution.  Then  the  latex  suspension  was 
pumped  through  the  column.  The  effluent  particle  concentration 
was  continuously  monitored  by  a  UV  spectrophotometer  equipped 
with  a  flow-through  cuvette  (Shimadzu  UV-160).  Linear  correla¬ 
tions  between  particle/tracer  concentration  and  UV  absorbance  at 
a  wavelength  of  220  nm  (for  particles)  and  305  nm  (for  tracer)  were 
determined  in  preliminary  experiments.  The  influent  was  switched 
to  the  buffer  solution  after  3-5  pore  volumes  of  the  particle  suspen¬ 
sion.  The  filtration  experiment  was  finally  stopped  after  flushing 
with  another  3-4  pore  volumes  of  buffer  solution.  Three  different 
flow  rates  were  employed  in  the  filtration  of  positively  charged  par¬ 
ticles:  6.2,  3.2,  and  1.5mL/min,  corresponding  to  Darcy  velocities 
of  17.6, 9.1, 4.3  m/d,  respectively.  For  the  negatively  charged  parti¬ 
cles,  we  tested  the  same  three  flow  rates.  The  differences  in  effluent 
concentration  for  these  three  flow  rates  were,  however,  relatively 
small  such  that  we  could  not  observe  how  flow  rate  affects  par¬ 
ticle  deposition.  To  better  resolve  this  dependency,  we  used  flow 
rates  of6.2, 1.5  and  0.37  mL/min,  corresponding  to  Darcy  velocities 
of  17.6,  4.3  and  1.1  m/d,  respectively,  for  the  negatively  charged 
particles.  Note  that  the  Darcy  velocities  used  in  the  experiments 
are  higher  than  those  typically  observed  for  groundwater  flow, 
but  they  can  be  representative  for  filtration  during  water  treat¬ 
ment.  During  each  filtration  experiment,  the  flow  rate  remained 
constant.  All  experiments  were  conducted  at  room  temperature, 
25  °C. 

2.5.  Micro-computed  tomography 

We  obtained  high-resolution  volumetric  images  of  the  filtration 
column  by  using  non-destructive  X-ray  micro-computed  tomogra¬ 
phy,  as  used  in  many  other  areas,  including  soil  science  |2]  and 
hydrology  [42].  We  used  a  microfocus  X-ray  computed  tomogra¬ 
phy  (MXCT)  system.  Our  MXCT  produces  polychromatic  X-rays 
with  an  acceleration  potential  ranging  from  10  to  225  keV  and 
an  electrical  current  ranging  from  0  to  3  mA.  During  a  CT  scan, 
X-rays  are  projected  from  a  cone-beam  onto  and  through  the  sam¬ 
ple,  which  is  located  between  the  X-ray  source  and  the  detector 
(a  charge-coupled  device  or  CCD).  The  sample  is  rotated  by  360‘ 
at  predetermined  angular  increments.  X-ray  attenuation  depends 
on  the  sample  density,  atomic  number,  and  column  diameter.  Pro¬ 
vided  a  sufficient  number  of  X-rays  pass  through  the  sample  to  the 
detector,  a  senes  of  two-dimensional  (2D)  horizontal  images  (or 
CT  “slices")  are  produced  using  the  Feldkamp  algorithm  [16]  from 
the  X-ray  attenuation  data,  which  was  collected  at  2400  angular 
increments  during  sample  rotation.  These  2D  images,  which  are 
accurately  aligned  in  the  vertical  dimension,  are  converted  into  a 
volumetric  image  during  post-processing.  The  cone-beam  captures 
the  full  diameter  of  the  sample  column,  but  because  the  height  of  a 
sample  column  is  often  longer  than  the  diameter  of  the  X-ray  beam, 
the  sample  is  translated  vertically  and  scanned  again.  The  process 
of  translating  the  sample  vertically  and  then  scanning  it  is  repeated 
until  the  desired  sample  length  is  imaged. 

3.  Theoretical  methods 

3. 1 .  Simulation  of  fluid  flow  by  the  LB  method 

To  simulate  filtration,  first  the  fluid  velocity  field  is  determined. 
At  the  pore  scale,  the  Navier-Stokes  equation, 

dv  -  --  1  -  -2-  ^ 

— — v  ■  V  i/  = - Vp-j-vV  v-tg 

3t  pf 


(2) 
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describes  fluid  flow,  where  v  is  the  pore-scale  velocity  field,  p  [Pa  ]  is 
pressure,  v  [m2/s]  is  the  kinematic  viscosity,  pf  [kg/m3]  is  the  fluid 
density,  and  g  [m/s2]  is  the  gravitational  acceleration  vector. 

The  LB  method  represents  the  fluid  particles  by  quasi  particles 
that  travel  along  a  typically  regular  numerical  lattice  and  collide  at 
the  lattice  nodes.  We  used  the  LB  model  described  in  [22]  for  solv¬ 
ing  the  Navier-Stokes  equation  (2)  numerically.  This  model  uses  a 
three-dimensional  cubic  lattice  with  15  velocity  vectors  (a  D3Q15 
model  with  6  velocities  along  the  lattice  axes,  8  along  the  space 
diagonals,  and  1  zero  vector  for  rest  particles).  The  LB  equation  for 
single-phase  flow  can  be  written  as 

/n(*i  +  Cn,  t -F  1 )  — /n(*i.  0  =  (3) 


As  we  do  not  account  for  anisotropic  diffusion,  the  diffusion  ten¬ 
sor  becomes  a  scalar  quantity  which  can  be  estimated  from  the 
Einstein-Stokes  equation, 

D=  kT  (8) 

3wpfvdp 

where  dp  [m]  is  the  diameter  of  the  suspended  spherical  particles. 

Our  LB  method  uses  the  same  numerical  lattice  for  particle 
transport  as  is  used  for  modeling  fluid  flow  (D3Q15).  The  LB  equa¬ 
tion  for  particle  transport  is 

C„(x,  +  e„,  t  +  1)  -  C„(x„  t)  =  -^0)(x„  t)  -  Cn(x„  I)]  (9) 

<-C 


where fn  is  the  distribution  function  for  the  fluid,  is  the  collision 
term,  and  n  =  0,  . . .  ,14  specifies  the  velocity  vector.  The  left  hand 
side  of  Eq.  (3)  describes  the  travel  (streaming)  step,  and  the  right 
hand  side  describes  the  collision  step.  We  use  a  single-relaxation¬ 
time  model  for  the  collision  term: 

i2n=^l^°)(Xi.t)-/„(xj.OI  (4) 

rf 

where  zy  is  the  relaxation  time  of  the  fluid  which  is  related  to  the 
fluid  viscosity  by 


r/  = 


1  v 

2  Axvq 


(5) 


where  Cn  [-]  is  the  distribution  function  for  the  suspended  particles, 
and  tc  [-]  is  the  relaxation  time  for  the  suspended  particles  which 
is  related  to  the  lattice  Peclet  number  Pec  by 

_L  °  *  1  00) 

Pec  Axvo  3  6 

If  the  external  force  in  Eq.  (7)  accounts  for  gravity  and  buoy¬ 
ancy,  then  the  nondimensional  advection-diffusion  equation  can 
be  approximated  by  the  LB  equation  (9)  with  an  equilibrium  distri¬ 
bution  function  given  by  [27] 

c'0,  =  /3„c[l  +  3e„(V  +  Vs)|  (ll) 


vq  -  Ax/At,  Ax  is  the  lattice  spacing,  At  is  the  time  step, 

r  a  9  a  2  3^2 

fi0)  =  ynp  l+3Ve„  +  j(l/  e„)  - 


(6) 


is  the  equilibrium  distribution  function  for  the  fluid,  yo  =  2/9, 

Yn“ 1/9  for  n-1 .  6,  yn - 1/72  for  7 .  14, 

is  the  nondimensional  fluid  density,  which  is  related  to  the 

nondimensional  pressure  by  P  =  p/3,  and  V  =  Yl]=ofi fo/p  IS 
nondimensional  fluid  velocity.  See  the  appendix  in  [22]  fora  deriva¬ 
tion  that  shows  that  the  LB  method  simulates  fluid  flow  fields  which 
satisfy  the  Navier-Stokes  equation. 

In  a  LB  simulation,  one  can  induce  flow  by  applying  either  a  body 
force  or  a  pressure  gradient  [8,10],  A  body  force  requires  periodic 
boundary  conditions  which  do  not  occur  in  the  experiment.  We 
generate  a  0.1%  pressure  (density)  difference  between  the  inlet  and 
the  outlet  to  induce  flow.  Fluid  flows  in  the  vertical  direction,  i.e. 
upward.  The  velocity  field  is  then  retrieved  from  the  simulations  at 
steady-state.  This  velocity  field  is  used  as  an  input  to  the  particle 
transport  model  described  next. 


32.  Simulation  of  panicle  transport  and  deposition  by  the  LB 
method 


An  advection-diffusion  equation  [11]  describes  particle  trans¬ 
port  at  the  pore  scale 

|£  +  v.(5c)=  V.(DpVc)-V.  (7) 

where  k  [J/K]  is  the  Boltzmann  constant,  T  [K]  is  the  absolute  tem¬ 
perature,  F  [N]  is  the  external  force  on  the  particles  (e.g.  gravity 
and  DLVO  forces),  and  Dp  [m2/s]  is  the  particle  diffusion  tensor. 
In  our  simulations,  gravity  is  the  only  external  force.  We  do  not 
account  for  the  lubrication  effect  [31  ]  and  anisotropic  diffusion  [  1  ], 
both  of  which  are  accounted  for  in  the  correlation  equation  by  [40] 
for  clean-bed  filtration  in  sphere  packings.  Even  though  this  latter 
correlation  was  developed  for  sphere  packings,  we  could  actually 
use  it  to  estimate  the  effects  of  lubrication  and  anisotropic  diffu¬ 
sion  on  particle  transport  in  our  packings  of  flattened  half  spheres. 
We  found  that  these  effects  should  affect  pD  only  marginally  [27 1. 


where  =  2/9,  m  1  /9  for  n  =  1 . 6,  =  1  /72  for  n  =  7,  ....  14, 

c  =  M  *s  nondimensional  particle  concentration, 


’  (Pp-Pf)zdj 
18  pfWo 


(12) 


is  the  nondimensional  particle  settling  velocity,  and  pp  [kg/m3]  is 
the  particle  density. 

Through  a  Chapman  and  Enskog  expansion,  we  can  show  that 
the  concentration  fields  obtained  from  the  LB  equation  (9)  with  the 
equilibrium  distribution  functions  given  by  Eq.  ( 1 1  )do  approximate 
the  ones  from  the  nondimensional  advection-diffusion  equation 
[20].  Moreover,  our  modeling  approach  [27]  accounts  for  all  of  the 
three  transport  mechanisms,  diffusion,  interception  and  gravita¬ 
tion,  which  according  to  [45],  contribute  to  particle  filtration. 


3.3.  Development  of  a  correlation  equation  for  rj0 


We  performed  numerical  column  breakthrough  experiments  in 
the  3D  pore  structure  of  the  volumetric  X-ray  CT  image  and  inferred 
r\D  from  the  steady-state  effluent  concentration  c0Uf.  We  first  sim¬ 
ulated  the  steady-state  pore-scale  velocity  field  in  the  scanned 
porous  medium.  Then,  we  used  this  velocity  field  as  an  input  to  our 
LB  model  for  particle  transport.  We  applied  the  Dinchlet  boundary 
condition  at  the  inlet  of  the  column,  c(x,y,z  =  0,  t)  =  cin,  where 
q„  is  the  concentration  of  the  suspended  particles  injected  into 
the  column.  At  the  outlet,  we  used  the  Neumann  boundary  condi¬ 
tion,  dc(x,y,z  =  L,  t)/dz  =  0.  The  perfect-sink  boundary  condition, 
c-0,  at  the  collector  surface  could  be  applied  to  obtain  a=  1  since 
dp  «  dc.  When  deriving  a  correlation  equation  for  t?0,  we  could 
assume  that  interception  and  gravitation  had  negligible  impact  on 
ho.  i.e.  r)D  ho.  because  particle  size  dp  was  very  small.  The  collector 
efficiency  ho  can  then  be  determined  from  [44,45] 


hD  =  - 


2  dc 
3(1  -e)L 


(13) 


Eq.  (13)  is  obtained  by  solving  Eq.  (1)  using  the  assumptions 
that  flow  is  one-dimensional  and  at  steady-state,  and  longitudinal 
hydrodynamic  dispersion,  which  occurs  at  the  Darcy  scale,  does  not 
significantly  affect  particle  transport  [45]. 
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Fig.  1.  CT  images  of  the  packing,  (a)  A  vertical  and  horizontal  slice  The  red  represents  the  solid  phase,  and  the  blue  represents  the  pore  space,  (b)  3D  view. 


For  nonspherical  grains  such  as  our  collector  particles,  the  sur¬ 
face  average  equivalent  sphere  diameter,  Dc  [m],  is  often  used  to 
specify  their  size: 


where  Vc  [m3]  is  the  grain  volume,  and  Sc  [m2]  is  the  surface  area 
of  the  grain.  We  can  determine  Vc  and  Sc  from  the  volumetric  CT 
images.  As  we  neglect  the  lubrication  effect  and  van  der  Waals 
forces  (we  do  not  include  Nvdw  in  our  correlation),  we  assume  the 
correlation  for  pD  to  be  a  function  of  the  porosity,/(f ),  the  intercep¬ 
tion  number,  Nr  =  dp/Dc,a nd  the  Peclet  number  NPe  -  C/Dc/D  where 
U  is  the  Darcy  velocity.  We  use  the  expression/^)  =  (1  -  ef  /e2  to 
model  the  porosity  dependence  of  r}D.  Such  a  porosity  dependence 
has  also  been  used  in  a  previous  study  [27)  on  the  collector  effi¬ 
ciency  of  packings  of  spheres.  The  concentration  limits,  couf  =  0  for 
e  *  0  and  c0ur  *  Qn  for  e  -*  1 ,  are  properly  captured  by  our  J[e).  Thus. 
tjD  can  be  written  as 

>,D  =  C,^^NgNjj3  (15) 

The  coefficients  C\ ,  C2,  and  c3  in  Eq.  ( 1 5)  can  be  determined  from 
a  multivariate  regression  analysis  on  the  r)o  values  obtained  from 
a  set  of  the  numerical  breakthrough  experiments  in  the  volumetric 
CT  image  of  the  filter  material. 

4.  Results  and  discussion 

4. 1.  Experimental  results  from  X-ray  CT 

We  imaged  a  middle  layer  of  the  entire  column.  The  MXCT  was 
operated  with  an  acceleration  potential  of  1 00  keV  and  an  electrical 
current  of  80  pA  The  sample  stage  was  set  65  cm  from  the  X-ray 
source.  This  setup  enabled  the  X-ray  beam  to  span  the  2.54  cm  in 
diameter  sample  and  allowed  X-rays  to  reach  the  detector  at  a  suf¬ 
ficient  intensity.  After  each  scan  of  the  sample,  the  sample  was 
translated  vertically  to  scan  a  desired  sample  height.  Given  this 
setup,  the  sample  was  imaged  in  the  x,  y,  and  z-directions  with  a 
resolution  of  0.05045  mm.  Given  the  number  of  voxels  per  particle 
diameter,  this  resolution  provides  an  accurate  representation  of  the 
PVDF  beads;  accurate  representation  of  similar  particles  has  been 
determined  to  occur  for  10-20  voxels  per  particle  diameter  with 
accuracy  increasing  as  the  number  of  voxels  per  particle  diameter 
increases  ] 39]. 

The  entire  image  of  the  middle  layer  consists  of  51 9  x  51 9  x  533 
voxels.  MXCT  images  are  16-bit  gray-scale  images.  A  good  image 
contrast  was  obtained  due  to  the  X-ray  attenuation  difference 


between  the  flattened  half  spheres  and  the  surrounding  air.  While 
the  gray-scale  histogram  of  the  MXCT  image  has  two  distinct  peaks, 
one  for  the  flattened  half  spheres  and  one  for  the  surrounding  pore 
space,  a  percentage  of  voxels  has  intermediary  values,  which  are 
difficult  to  assign  to  a  given  population  using  a  simple  thresholding 
technique.  Therefore,  an  indicator  kriging  algorithm  was  employed 
to  convert  the  entire  MXCT  image  into  a  1-bit  binary  image.  This 
algorithm  is  used  to  assign  voxels  that  clearly  belong  to  the  flat¬ 
tened  half  spheres  a  value  of  1  and  to  assign  voxels  that  clearly 
belong  to  the  pore  space  a  value  of  0.  The  intermediary  values  are 
assigned  a  value  of  1  or  0  depending  upon  the  gray-scale  values  of 
the  neighboring  voxels  (28).  Porosity  is  then  readily  determined  by 
taking  the  sum  of  the  0’s  divided  by  the  total  number  of  voxels  in 
the  image.  Fig.  1  shows  that  the  resulting  image  provides  a  realistic 
geometry  that  resolves  the  orientation  of  the  PVDF  beads  and  also 
the  contact  points  between  the  beads. 

Through  image  analysis,  we  determined  that  the  middle  layer 
had  a  porosity  of  38%.  This  porosity  is  lower  than  the  average  col¬ 
umn  porosity  (43%),  indicating  that  the  beads  were  not  uniformly 
packed.  Ideally,  we  should  have  scanned  the  entire  column  to  iden¬ 
tify  a  region,  in  which  the  porosity  was  about  equal  to  the  average 
column  porosity  of  43%,  and  then  used  this  region  in  the  filtration 
simulation.  It  was,  however,  unrealistic  to  scan  the  entire  column 
due  to  the  cost  of  the  MXCT  measurements.  In  order  to  account  for 
the  difference  in  porosity,  we  used  the  porosity  of  the  middle  layer 
when  determining  the  parameters  of  the  correlation  equation  (15) 
from  the  LB  simulations,  because  this  porosity  corresponds  to  the 
simulation  domain.  When  predicting  po  of  the  entire  column,  how¬ 
ever,  we  used  the  porosity  of  the  entire  column  in  our  correlation 
equation. 

4.2.  Results  for  deposition  experiments 

Using  the  strategy  explained  in  [14],  we  experimentally  deter¬ 
mined  conditions  for  which  a-1.  The  idea  is  that  any  particle 
colliding  with  a  collector  will  attach  to  it,  if  the  collectors  and  the 
particles  are  oppositely  charged.  Such  experiments  are  performed 
at  a  relatively  high  ionic  strength  in  order  to  avoid  ot>  1,  that  is, 
deposition  enhancement  due  to  attractive  van  der  Waals  forces  [9]. 
To  prevent  particle  aggregation,  excessively  high  ionic  strengths 
(>5  mM)  are  avoided. 

We  conducted  particle  deposition  experiments  with  the  posi¬ 
tively  charged  Latex  A  particles.  Since  pristine  PVDF  collectors  are 
possibly  only  slightly  negatively  charged  due  to  the  lack  of  ion- 
izable  function  groups  on  PVDF  structures,  coating  the  collectors 
with  the  surfactant  SDS,  which  possesses  a  negatively  charged  sul¬ 
fate  group,  should  increase  the  zeta  potential  and  the  surface  charge 
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Fig.  2.  Particle  breakthrough  curves  for  Latex  A  particles:  (a)  PVDF  beads  either 
coated  or  not  coated  with  surfactant(SDS).  (b)  PVDF  beads  are  coated  withSDS;  use 
of  three  different  flow  rates. 

density  of  the  collectors  |24|.Todetermine  the  conditions  for  which 
a  =  1,  we  obtained  breakthrough  curves  with  the  collectors  either 
being  coated  with  SDS  or  not.  We  used  a  CaCh  solution  with  an 
ionic  strength  of  1  mM  to  prevent  a>  1.  The  breakthrough  curves 
in  Fig.  2(a)  show  that  a  is  higher  for  the  SDS  coated  collectors.  Thus, 
we  concluded  that  the  SDS  coated  collectors  yielded  a  =  1 .  Using  the 
SDS  coated  collectors,  a  set  of  triplicate  deposition  experiments  was 
performed  for  three  different  flow  rates.  Fig.  2(b)  illustrates  the  flow 
rate  dependence  of  the  breakthrough  curves. 

4  3.  Modeling  results 


Fig  4.  Numerical  breakthrough  curve  for  the  simulation  illustrated  in  Fig  3. 


[  1 8  ]  to  determine  the  average  surface  area  of  a  collector,  Sc,  from  the 
CT  image.  Likewise,  the  collector  volume,  Vc,  was  calculated  from 
the  segmented  CT  image.  Using  Eq.  (14),  we  obtained  Dc-  3.6  mm. 

Because  of  the  computational  expense  filtration  was  simulated 
in  subdomains  of  the  volumetric  CT  image.  Due  to  the  random¬ 
ness  of  the  pore  structure,  r)D  differs  for  different  subdomains  of 
the  porous  medium,  even  if  porosity  e  is  the  same,  as  we  also 
found  in  a  study  on  filtration  in  sphere  packings  (27).  The  mean 
value  of  t)D  is  expected  to  remain  approximately  constant,  while 
the  variation  in  r\D  is  expected  to  decrease,  as  the  domain  size 
increases.  The  domain  size  has  reached  a  representative  elementary 
volume  (REV),  once  the  variation  becomes  negligible.  To  account 
for  the  possibility  that  the  volumes  of  the  subdomains  are  smaller 
than  a  REV,  we  identified  four  different  subdomains,  all  of  which 
have  about  the  same  porosity,  e  =  0.38  ±0.01.  The  subdomains 
consist  of  250  x  250  x  200,  250  x  250  x  200,  200  x  200  x  300,  and 
200  x  200  x  400  voxels  and  contain  between  80  and  105  flattened 
half-spherical  collectors.  The  subdomains  have  slightly  different 
volumes  in  order  to  ensure  that  their  four  e  values  fall  into  the 
prescribed  e  interval. 

In  these  four  subdomains,  we  performed  a  set  of  numerical 
breakthrough  experiments  for  different  values  of  U  and  dp.  Fig.  3 
illustrates  how  the  particle  concentration  develops  during  such  a 
numerical  experiment.  The  breakthrough  curve  is  shown  in  Fig.  4. 
For  each  value  of  U  and  dp,  we  obtained  four  different  tjD  values, 
from  which  we  calculated  an  average  i)d  value.  Those  average  val¬ 
ues  were  then  used  to  develop  a  correlation  equation  for  i]D.  We 
determined  the  coefficients  in  Eq.  (15)  through  regression  analysis 
and  obtained 


To  determine  the  coefficients  in  Eq.  (15)  we  first  inferred  Dc  from 
the  volumetric  CT  image.  We  used  a  pore-space  analysis  software 


rjD  =(8.34  ±0.75) 


a-£)3 

£2 


jw  0.S3±0.01S 
nPe 


M  0.03±0.017S 


(16) 


Fig  3.  Particle  concentration  at  selected  times  t  in  a  packing  of  flattened  half  spheres.  The  color  represents  the  particle  concentration  and  the  black  represents  the  solid 
phase.  Note  that  fluid  flows  from  z- Otoz-  1.01  cm  while  the  simulation  results  are  displayed  oppositely  to  give  a  better  visualization.  Simulation  parameters:  U- 0.025  [xm/s. 
£-0.38.  dp  -20  nrn.and  L- 1.01  cm. 
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Fig.  5.  Comparison  of  the  collector  efficiency  from  the  numerical  simulations  with 
the  predictions  from  the  correlation  for  the  four  different  subdomains  and  for  dif¬ 
ferent  values  of  U  and  dp.  Simulation  parameters:  Dc  -  3.6  mm,  e  -  0.43,  L -  1.01  mm. 
and  Ax«S0.48pm. 


Table  2 

Summary  of  parameters  used  in  numerical  simulations. 


Parameter 

Range 

0.02-1  pm 

e 

0.38 

U 

0.02S-1  pm/s 

T 

298  K 

k 

1.38  x  10‘23J/K 

Dc 

3.6  mm 

Pp 

1.0SSg/cm3 

Fig.  5  compares  the  collector  efficiency  inferred  via  Eq.  ( 1 3 )  from 
individual  numerical  breakthrough  experiments,  rjsimt  with  the  pre¬ 
dictions  based  on  the  correlation  equation  (16),  r)COrr .  The  good 
agreement  between  r}sjm  and  qcorr  demonstrates  the  appropriate¬ 
ness  of  the  functional  form  of  our  correlation  for  rjD.  A  correlation 
coefficient  of  R  =  0.98  illustrates  the  good  fit.  Table  2  summarizes 
the  parameters  used  in  the  simulations. 

Note  that  the  velocities  we  used  in  our  simulation  vary  from 
0.002  to  0.1  m/d  and  are  smaller  than  the  ones  we  used  in  the  exper¬ 
iments.  This  is  due  to  numerical  instability  which  would  occur  for 
high  Peclet  numbers  or  high  flow  rates.  Due  to  the  relatively  small 
Peclet  numbers  used  in  the  simulations,  the  simulated  nondimen- 
sional  breakthrough  concentrations  are  small  too  (see  Fig.  4).  Even 
though  the  correlation  equation  that  we  developed  for  r}D  is  based 
on  simulations  with  small  Peclet  numbers,  the  correlation  equation 
can  be  used  to  predict  filtration  experiments  performed  at  higher 
Peclet  numbers  as  demonstrated  in  a  similar  study  on  filtration  in 
sphere  packings  [27].  As  in  [27]  we  will  use  statistical  methods  to 
quantify  the  extrapolation  error. 

We  can  also  comment  on  the  REV  for  filtration  in  a  packing  of 
flattened  half  spheres,  even  though  this  study  has  not  been  specif¬ 
ically  designed  to  determine  this  REV.  To  that  end,  we  would  have 
needed  to  simulate  filtration  for  a  wide  range  of  domains  sizes, 
as  in  LB  studies  of  permeability  [5,46,29]  or  in  pore-morphology 


Comparison  of  the  Experiments  with  the  Correlation 


Fig.  6.  Comparison  of  the  collector  efficiency  from  the  correlation  and  the  experi¬ 
ments  for  Latex  A  particles  for  three  different  flow  rates:  6.2  mL/min  ( v ).  3.2  mL/min 
CO),  and  l.S  mL/min  (A). 

based  modeling  studies  of  primary  drainage  [21 1.  We  can,  however, 
examine  the  statistics  of  r]D  for  the  simulations  in  the  four  statisti¬ 
cally  identical  subdomains.  A  possible  definition  of  REV  relates  to 
the  domain  volume,  at  or  above  which  the  relative  standard  devi¬ 
ation  of  the  variable  of  interest  is  smaller  than  a  user-specified 
percentage,  for  example,  5%  in  the  permeability  study  by  [29]. 
Table  3  shows  that  the  relative  standard  deviation  of  r\D  is  a  func¬ 
tion  of  both  Darcy  velocity  U  and  particle  size  dp.  The  relative 
standard  deviation  increases  with  decreasing  dp  and  increasing 
U  values.  Our  simulation  domain  would  constitute  a  REV  for  all 
dp  and  U  values  presented,  if  one  accepted  a  10%  standard  devia¬ 
tion.  If  one  accepted  a  5%  standard  deviation,  then  our  simulation 
domain  would  be  a  REV  only  for  dp  above  200  nm  and  U  below 
2.5  x  10  8  m/s.  Hence,  depending  on  the  accepted  standard  devi¬ 
ation,  our  simulation  domains  may  constitute  a  REV  or  not.  We 
accounted  for  the  latter  possibility  by  having  performed  filtration 
simulations  in  four  different  subdomains  of  the  imaged  packing  of 
flattened  half  spheres  that  all  have  the  same  porosity  f.  Simula¬ 
tions  in  bigger  domains  would  be  desirable  and  are  expected  to 
yield  smaller  error  bounds  for  the  coefficients  in  our  correlation 
equation  (16). 

Note  that  our  correlation  equation  given  by  Eq.  (16)  accounts 
explicitly  for  the  porosity  e.  However,  the  correlation  equation  is 
based  on  simulations  in  packings  with  a  single  porosity  of  38%.  We 
have  neither  experimental  nor  numerical  data  to  support  the  poros¬ 
ity  dependence  posited  in  Eq.  (16).  Additional  investigations  are 
necessary  to  test  the  porosity  term  in  Eq.  (16).  We  did  not  perform 
these  investigations  due  to  the  cost  of  obtaining  CT  scans. 

4.4.  Comparison  of  rjD  between  the  correlation  equation  and  the 
experiments 

Fig.  6  compares  the  collector  efficiency  measured  in  the  exper¬ 
iments,  rjexp,  with  the  prediction  calculated  from  our  correlation 
equation  (16),  rjcon-  The  rjexp  values  reflect  the  mean  values  from 
the  triplicate  deposition  column  experiments  that  we  performed 
for  the  three  different  flow  rates.  In  order  to  evaluate  the  pre- 


Table  3 

Statistics  of  the  r}D  values  simulated  in  statistically  identical  packings  of  flattened  half  spheres. 


dp  (nm) 

Relative  standard  deviation  of  r}D 
(constant  U- 2.S  x  10"8  m/s) 

U(mfs) 

Relative  standard  deviation  of  tiD 
(constant  dp-l00nm) 

20 

8.42% 

2.47x10® 

4.9  S% 

SO 

7.86% 

6.17x10  8 

7.S1% 

100 

S.72% 

1.23  x  10-7 

8.72% 

200 

1.4S% 

2.47  x  10  7 

9.93% 
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Table  4 

Removal  efficiency  for  three  different  nanoparticle  sizes  and  three  different  flow 
rates  in  deposition  experiments. 


Particle  label 

Flow  rate(mL/min) 

Porosity.  f(%) 

1  —CoutlCin  (%) 

Of  [-1 

Latex  B 

0.35 

43 

79 

0.24 

Latex  B 

1.5 

43 

63 

0.32 

Latex  B 

6.15 

43 

46 

0.42 

Latex  C 

0.39 

43 

71 

0.2B 

Latex  C 

1.53 

43 

55 

0.3B 

Latex  C 

6.15 

43 

40 

0.51 

Latex  D 

0.37 

43 

70 

0.25 

Latex  D 

1.5 

43 

45 

0.27 

Latex  D 

6.1 

43 

34 

0.39 

dictive  capability  of  Eq.  (16).  we  quantified  the  uncertainties  of 
both  the  prediction  from  the  correlation  equation  and  the  mea¬ 
sured  values.  The  smaller  of  the  two  overlapping  vertical  error  bars 
represents  the  95%  confidence  interval  of  the  prediction  of  the  cor¬ 
relation  equation.  In  statistical  terms,  this  is  the  confidence  interval 
for  the  mean  response  |41,  Ch.  12].  If  the  experimental  measure¬ 
ment  was  without  error,  then  this  error  bar  should  intersect  with 
the  one-to-one  line  for  the  correlation  to  explain  the  experimental 
data.  This  is  the  case  for  all  the  three  data  points.  The  mean  experi¬ 
mental  riexp  values  have,  however,  some  error,  which  we  quantified 
by  the  95%  confidence  interval  of  the  mean.  Hence,  the  prediction 
of  the  correlation  equation  may  deviate  from  the  one-to-one  line 
by  this  confidence  interval.  In  other  words,  the  sum  of  the  errors 
from  the  correlation  and  the  experimental  measurement  should 
intersect  with  the  one-to-one  line.  The  bigger  of  the  two  overlap¬ 
ping  vertical  error  bars  represents  this  combined  error.  As  these 
combined  errors,  and  even  the  confidence  intervals  of  the  mean 
response  alone,  overlap  with  the  one-to-one  line,  we  can  conclude 
that  the  correlation  agrees  well  with  the  results  of  the  column 
experiments.  Correlation  equations  developed  for  spherical  collec¬ 
tors  [45,40,33]  do  not  and  cannot  expected  to  agree  well  with  our 
experimental  measurements  (companson  not  presented)  because 
of  the  difference  between  spherical  and  our  about  half-spherical 
collectors. 

Our  comparison  between  the  experimental  data  and  our  cor¬ 
relation  equation  is  hampered  by  the  fact  that  our  column  was 
not  packed  uniformly  as  exemplified  by  a  spatially  variable  poros¬ 
ity.  Hence,  the  measured  collector  efficiency,  ^exp.  should  depend 
not  only  on  porosity  but  also  on  the  porosity  variation.  We  could 
not  quantify  the  effect  that  a  non-uniform  porosity  has  on  t]0,  but 
we  expect  the  dependence  of  r)o  on  the  porosity  variation  to  be  a 
correction  to  the  dependence  on  e ,  NPe,  and  NR.  Additional  inves¬ 
tigations,  which  would  involve  scanning  entire  filtration  columns 
and  testing  them  for  uniform  porosity,  are  necessary  to  resolve  this 
issue. 

We  also  examined  whether  our  assumption  of  clean-bed  filtra¬ 
tion  was  justified.  To  that  end,  we  assumed  that  all  particles  are 
deposited  onto  the  collector  surfaces  during  the  entire  filtration 
event  and  estimated  the  volume  of  particle  deposited.  Using  a  mass 
balance,  we  calculated  a  reduction  in  porosity  by  0.00019%  when 
flushing  a  column  with  one  pore  volume  of  nanoparticle-laden 
solution.  In  the  experiments,  the  columns  were  flushed  with  4-6 
pore  volumes.  Hence,  the  porosity  was  reduced  at  most  by  about 
0.001%.  This  small  number  justifies  the  assumption  of  clean-bed 
filtration. 

4.5.  Application  of  the  correlation  to  predict  a 

The  negatively  charged  Latex  B,  C,  and  D  particles  (see  Table  1 ) 
used  together  with  the  negatively  charged  collectors  (without  SDS 
coating)yield  unfavorable  chemical  conditions  with  respect  to  par¬ 
ticle  deposition  (significantly  lowera  values).  For  these  unfavorable 


conditions,  the  collision  efficiency  a  was  calculated  as  follows  [45]: 


1  2  dc 
lo  3  (1  -e)L 


(17) 


where  Eq.  (16)  was  employed  for  r]o .  Table  4  summarizes  a  and 
the  removal  efficiency  1  -coutlcin  for  the  three  particle  sizes  and 
three  different  flow  rates.  As  predicted  by  colloid  filtration  theory 
[45],  the  removal  efficiency  increases  as  the  flow  rate  and  particle 
size  decrease,  because  in  this  nanoparticle  experiments,  particle 
deposition  is  governed  by  diffusion  while  gravitational  settlement 
is  negligible. 


Acknowledgments 

We  thank  Professor  Charles  R.  O’Melia  for  discussion  and  direc¬ 
tion,  Professor  Hugh  Ellis  for  computing  resources  at  Johns  Hopkins 
University,  Dr.  Roland  Glantz  for  his  help  on  the  calculation  of 
surface  area,  Atofina/Arkema  chemical  company  for  donating  the 
Kynar  samples  used  in  the  experiments,  and  Kevin  Wilson  for 
assistance  with  image  processing.  We  also  thank  the  ONR  Coastal 
Geosciences  Program  for  funding  for  the  CT  work.  MH  was  sup¬ 
ported  by  NSF  Grant  EAR-0911425. 

References 

[1]  Z.Adamczyk.T.  Dabros.J.Czarnecki.T.G.M.Vandeven,  Particle  transfer  to  solid- 
surfaces.  Advances  in  Colloid  and  Interface  Science  19  (1983)  1B3-252. 

|2]  S.H.  Anderson.  R.L  Peyton,  C.J.  Gantzer,  Evaluation  of  constructed  and  natural 
soil  macropores  using  X-ray  computed  tomography.  Geoderma  46  ( 1  -3)  (1 990) 
1 3-29. 

[31  R.C.  Bales.  S.M.  Li,  K.M.  Maguire.  M.T.  Yahya.  C.P.  Gerba,  Ms-2  and  poliovirus 
transport  in  porous  media  -  hydrophobic  effects  and  chemical  perturbations. 
Water  Resources  Research  29  (4)  (1993)  957-963. 

[4[  T.  Baumann.C.J  Werth.  Visualization  and  modeling  of  polystyrol  colloid  trans¬ 
port  in  a  silicon  micromodel.  Vadose  Zone  Journal  3  (2)  (2004)434-443. 

(5]  W  J.Bosl.J.Dvorkin,  A.  Nur,  A  study  of  porosity  and  permeability  using  a  Lattice- 
Boltzmann  simulation.  Geophysical  Research  Letters  25  (9) (1 998)  1475-1478. 

[6[  V.N.  Burganos,  E.D.  Skouras,  CA.  Paraskeva.  A.C.  Payatakes,  Simulation  of  the 
dynamics  of  depth  filtration  of  non-Brownian  particles.  AlChE  Journal  47  (4) 
(2001)BB0-894. 

{7]  C.  Chen,  A.L  Packman.  J.F.Gaillard,  Pore-scale  analysis  of  permeability  reduction 
resulting  from  colloid  deposition.  Geophysical  Research  Letters  35  (7)  (200B). 
art.  no.  L07404. 

[B]  H,  Chen.  S.  Chen.  W.H.  Matthaeus.  Recovery  of  the  Navier-Stokes  equa¬ 
tions  using  a  lattice-gas  Boltzmann  methods,  Physical  Review  A  45  (1992) 
5339-5342. 

|9[  K.L  Chen.  S.E.  Mylon.  M.  Elimelech,  Aggregation  kinetics  of  alginate-coated 
hematite  nanoparticles  in  monovalent  and  divalent  electrolytes.  Environmen¬ 
tal  Science  8i  Technology  40(2006)  1516-1523. 

1 101  S.  Chen,  C.D  Doolen.  Lattice-Boltzmann  method  for  fluid  flows.  Annual  Review 
of  Fluid  Mechanics  30  ( 199B)  329-364. 

(Ill  M.  Elimelech,  Particle  deposition  on  ideal  collectors  from  dilute  flowing  sus¬ 
pensions:  mathematical  formulation,  numerical  solution,  and  simulations. 
Separations  Technology  4  (1994)  186-212. 

[12]  M.  Elimelech,  C.R.  O’Melia,  Kinetics  of  deposition  of  colloidal  particles  in  porous 
media.  Environmental  Science  &  Technology  24  (10)  (1990)  152B-1536. 

[13[  M.  Elimelech,  C.R.  O’Melia.  Effect  of  particle  size  on  collision  efficiency  in  the 
deposition  of  Brownian  particles  with  electrostatic  energy  barriers.  Langmuir 
6  (6)  (1990)  1153-1 163. 

[14]  M.  Elimelech.  J.  Gregory,  X.  Jia,  R.  Williams,  Modeling  of  particle  deposition 
onto  ideal  collectors,  in:  RA  Williams  (Ed.),  Particle  Deposition  and  Aggrega¬ 
tion,  Butterworth-Heinemann  Ltd.,  Linacre  House, Jordan  Hill,  Oxford  0X2  BDP, 
1995.  pp.  113-154. 

[  1 5  j  D.W.  Elliott.  W.X.  Zhang,  Field  assessment  of  nanoscale  biometallic  particles  for 
groundwater  treatment.  Environmental  Science  and  Technology  35  (24)  (2001 ) 
4922-4926. 

[16]  LA  Feldkamp,  LC.  Davis.  J.W.  Kress.  Practical  cone-beam  algorithm.  Journal  of 
the  Optical  Society  of  America  1  (6)  (1984)  612-619. 

117]  T.R.  Ginn.  B.D.  Wood.  K.E.  Nelson,  T.D.  Scheibe.  E.M.  Murphy.  TP.  Clement. 
Processes  in  microbial  transport  in  the  natural  subsurface.  Advances  in  Water 
Resources  25  (2002)  1017-1042. 

[IB]  R.  Glantz.  M.  Hilpert,  Tight  dual  models  of  pore  spaces.  Advances  in  Water 
Resources  31  (200B)  787-806. 

[19]  J.  Happel,  Viscous  flow  in  multiparticle  systems-slow  motion  of  fluids  relative 
to  beds  of  spherical  particles.  AlChE  Journal  4  (195B)  197-201 . 

[20]  M.  Hilpert,  Lattice-Boltzmann  model  for  bacterial  chemotaxis.  Journal  of  Math¬ 
ematical  Biology  51  (2005)  302-332. 

121  ]  M.  Hilpert.  C.T.  Miller.  Pore-morphology- based  simulation  of  drainage  in  totally 
wetting  porous  media.  Advances  in  Water  Resources  24  (3-4)  (2001 )  243-255. 


I 


W.  Long  er  al.  /  Colloids  ond  Surfaces  A:  Physicochem.  Eng.  Aspects  3 58  (2010)  163-171  171 


[22 [  S.  Hou.  Q  Zou.  S.  Chen.  G.  Doolen.  A.C.  Cogley.  Simulation  of  cavity  flow  by 
the  Lattice- Boltzmann  method.  Journal  of  Computational  Physics  118  (1995) 
329-347. 

[23]  Q.J.  Kang.  P.C.  Lichtner.  D.X.  Zhang.  Lattice-Boltzmann  pore-scale  model  for 
multicomponent  reactive  transport  in  porous  media.  Journal  of  Geophysical 
Research-Solid  Earth  111  (B5) (2006). 

[24]  K.J.  Kim.  A.G.  Fane.  M.  Nystrom,  A.  Pihlajamaki.  W,R.  Bowen.  H.  Mukhtar.  Eval¬ 
uation  of  electroosmosis  and  streaming  potential  for  measurement  of  electric 
charges  of  polymeric  membranes.  Journal  of  Membrane  Science  1 16  (2)  (1996) 
149-159. 

[25  [  X.Q.  Li.  D.W  Elliott.  W.X  Zhang.  Zero-valent  iron  nanopartides  for  abatement  of 
environmental  pollutants:  materials  and  engineering  aspects.  Critical  Reviews 
in  Solid  State  and  Materials  Sciences  31  (4)  (2006)  111-122. 

[26[  C.L  Lin.  J.D.  Miller.  Pore  structure  analysis  of  particle  beds  for  fluid  trans¬ 
port  simulation  during  filtration.  International  Journal  of  Mineral  Processing 
73  (2-4)  (2004)  2B 1-294. 

[27  [  W.  Long.  M  Hilpert.  A  correlation  for  the  collector  efficiency  of  Brownian  parti¬ 
cles  in  dean-bed  filtration  in  sphere  packings  by  a  Lattice-Boltzmann  method. 
Environmental  Science  and  Technology  43  (2009)4419-4424. 

[2B]  W.  Oh.  W.B.  Lindquist,  Image  thresholding  by  indicator  kriging.  IEEE  Transac¬ 
tions  on  Pattern  Analysis  and  Machine  Intelligence  21  (1999)  590-602. 

[29]  C.  Pan.  M.  Hilpert,  CT.  Miller.  Pore-scale  modeling  of  saturated  permeabilities 
in  random  sphere  packings.  Physical  Review  E  64  (6)  (2001 )  066702. 

[30]  C.  Pan.  M.  Hilpert.  C.T.  Miller.  Lattice-Boltzmann  simulation  of  two-phase  flow 
in  porous  media.  Water  Resources  Research  40  (2004)  W01501. 

[31  [  D.C.  Prieve.  M.M.J.  Lin,  Adsorption  of  Brownian  hydrosols  onto  a  rotating-disk 
aided  by  a  uniform  applied  force.  Journal  of  Colloid  and  Interface  Science  76 
(1980)32-47. 

(32 [  R.  Przekop.  A.  Moskal.  L  Gradon.  Lattice-Boltzmann  approach  for  description 
of  the  structure  of  deposited  particulate  matter  in  fibrous  filters.  Journal  of 
Aerosol  Science  34  (2)  (2003)  133-147. 

[33[  R.  Rajagopalan,  C.  Tien.  Trajectory  analysis  of  deep-bed  filtration  with  the 
sphere- in-cell  porous  media  model.  AIChEJournal  22  (1976)  523-533. 

[ 34 [  LLC.  Rehmann,  C.  Welty,  Stochastic  analysis  of  virus  transport  in  aquifers. 
Water  Resources  Research  35  (7)  (1999)  19B7-2006. 


[35]  J.E.  Saiers.  J.N,  Ryan.  Colloid  deposition  on  non-ideal  porous  media:  the  influ¬ 
ence  of  collector  shape  and  roughness  on  the  single-collector  efficiency. 
Geophysical  Research  Letters  32  (21 )  (2005),  art.  no.  L21406. 

[36[  N  Savage.  MS.  Diallo.  Nanomaterials  and  water  purification:  opportunities  and 
challenges.  Journal  of  Nanoparticle  Research  7  (4-5)(2005)  331-342. 

[37[  Y.  Sim,  C.V  Chrysikopoulos.  Three-dimensional  analytical  models  for  virus 
transport  in  saturated  porous  media.  Transport  in  Porous  Media  30  (1)(199B) 
B7-112. 

[3B]  E.D,  Skouras.  V.N.  Burganos.  C.A.  Paraskeva,  A.C.  Payatakes.  Simulation  of  the 
dynamic  behavior  of  horizontal  granular  filters.  Separation  and  Purification 
Technology  56(3)  (2007)  325-339. 

[39]  K.  E  Thompson.  CS.  Willson.  W.  L  Zhang.  Quantitative  computer  reconstruction 
of  particulate  materials  from  microtomography  images.  Powder  Technology 
163  (3)  (2006)  169-1B2. 

[40[  N.  Tufenkji,  M.  Elimelech,  Correlation  equation  for  predicting  single-collector 
efficiency  in  physicochemical  filtration  in  saturated  porous  media.  Environ¬ 
mental  Science  &  Technology  3B  (2004)  529-536. 

[41  ]  R.E.  Walpole.  R.H.  Myers.  Probability  and  Statistics  for  Engineers  and  Scientists. 
5th  ed.,  Macmillan  Publishing  Company.  New  York.  1993. 

[42 [  D.  Wildenschild.  J.W.  Hopmans,  C.M.P.  Vaz,  M.L  Rivers.  D.  Rikard.  B.S.B.  Chris¬ 
tensen.  Using  X-ray  computed  tomography  in  hydrology:  systems,  resolutions, 
and  limitations.  Journal  of  Hydrology  267  (3-4)  (2002 )  285-297. 

[43]  D.  Wolf-Gladrow.  A  Lattice-Boltzmann  equation  for  diffusion.  Journal  of  Statis¬ 
tical  Physics  79  (5/6)  (1995)  1023-1 032. 

[44]  K.M.  Yao.  Influence  of  suspended  particle  size  on  the  transport  aspect  of  water 
filtration.  PhD  thesis.  University  of  North  Carolina  at  Chapel  Hill,  196B. 

[45[  K.M.  Yao.  M.M.  Habibian,  C.R.  O’Melia.  Water  and  waste  water 
filtration-concepts  and  applications.  Environmental  Science  8?  Technology  5 
(1971)  1105. 

[46]  D  Zhang.  R.  Zhang.  S.  Chen.  W  E.  Soli.  Pore  scale  study  of  flow  in  porous  media: 
scale  dependency.  REV.  and  statistical  REV.  Geophysical  Research  Letters  27(B) 
(2000)  1195-1 19B. 

[47[  W.X.  Zhang.  Nanoscale  iron  particles  for  environmental  remediation:  an 
overview.  Journal  of  Nanoparticle  Research  5  (3-4)  (2003)  323-332 


